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ABSTRACT 



I review the properties of three X-ray accreting milHsec- 
ond pulsars observed with INTEGRAL. Out of seven re- 
cently discovered accretion-powered pulsars (one discov- 
ered by INTEGRAL), three were observed with the IN- 
TEGRAL satellite up to 300 keV. Detailed timing and 
spectral results will be presented, including data ob- 
tained during the most recent outburst of the pulsar HETE 
J 1900. 1-2455. Accreting X-ray millisecond pulsars are 
key systems to understand the spin and accretion history 
of neutron stars. They are also a good laboratory in which 
to study the source spectra, pulse profile, and phase shift 
between X-ray pulses in different energy ranges which 
give additional information of the X-ray production pro- 
cesses and emission environment. 

Key words: pulsars: individual (XTE J 1807-294, IGR 
J00291H-5934, HETE J1900.1-2455) - starts: neutron - 
X-ray: binaries. 



1. INTRODUCTION 



Many low-mass X-ray binaries (LMXB) consist of a neu- 
tron star accreting from a low mass companion star (< 
1 Mq). Accreting matter may spin up the neutron star 
(NS); therefore, one of the possible endpoints of the evo- 
lution of a low-mass X-ray binary is expected to be a mil- 
lisecond pulsar, i.e. a rapidly spinning NS with a rather 
weak, ^ 10* Gauss, surface magnetic field. Although ev- 
idence for rapidly spinning NS in LMXBs was inferred 
from the burst oscillations that were seen during type I 
X-ray bursts in several systems [see 22:, for a review], the 
detection of millisecond pulsations in persistent emission 
remained elusive for many years until the discovery of 
the first accreting millisecond pulsar by Wijnands & van 
der Klis [27]. Since that time, a total of seven accreting 
MSPs have been detected. 

All of the accreting MSPs are X-ray transients; they 
spend most of the time in a quiescent phase, with X-ray 
luminosities of order of 10^^ — 10'^'^ erg s^^. They some- 
times show X-ray outbursts reaching X-ray luminosities 
of 10^^ — 10"^^ erg s^^, during which coherent pulsations 



are observed with frequencies in the range between 180 
and 435 Hz [see reviews by 28, 21]. This frequency is 
interpreted as the NS rotation frequency given by a hot 
spot (or spots) in an atmospheric layer of the rotating NS 
10] • MSPs represent a new class of objects connecting 
accretion powered X-ray pulsars with rotation powered 
millisecond radio pulsars. 

MSP energy spectra are successfully fitted by a two- 
component model consisting of a multicolor blackbody 
soft X-ray emission and a Comptonized spectrum for the 
hard X-ray emission. The soft thermal component could 
be associated with the radiation from the accretion disc 
and/or the heated NS surface around the shock [see e.g. 
[T6l I2TI1 . The hard emission is likely to be produced by 
thermal Comptonization in the hot accretion shock on the 
NS surface lisl I20I1 with the seed photons coming from 
the stellar surface. The observed hard spectra are sim- 
ilar to the spectra observed from atoll sources in their 
hard, low-luminosity state [1]. A representative spec- 
trum of the three observed MSP (XTE J1807-294 [6], 
IGR J0029 1 H-5934 L7J , and HETE J 1 900. 1 -2455 [8J ) with 
INTEGRAL is shown in Figure 1, 2 and 3. 

Recent years have shown INTEGRAL gave additional 
contributions to the study of MSPs. Our timing analy- 
sis of IGR J00291 H-5934 showed for the first time that 
INTEGRAL/ISGRI is capable of detecting the pulse pro- 
file of a 1.67 ms pulsar up to 150 keV. This allowed us to 
study the pulsed fraction up to this high energy. We con- 
firmed for the first time the increase of the pulsed fraction 
with energy in an accretion powered MSP, explained by 
the action of the Doppler effect on the exponentially cut- 
off Comptonization spectrum from the hot spot or from a 
componization model where we account the pulsed frac- 
tion to be produced in a corona cloud. We also measured 
soft time lags with a complex energy dependence. Simi- 
larly to SAX J 1808.4-3658 and XTE J 175 1-350, the time 
lags increase rapidly with energy until 10 keV. However, 
in IGR J00291 H-5934, the time lags first increase and then 
decrease slightly, saturating above 15 keV, and possibly 
reaching zero around 50 keV. This is a intriguing result, 
and if confirmed or observed in another object could be 
a serious challenge to any time-lag model. Additionally, 
with RXTE we measured an increase of the NS spin pe- 
riod during accretion for the first time. This provided a 
strong confirmation of the theory of 'recycled' pulsars in 



which the old neutron stars in LMXBs become milHsec- 
ond radio pulsars through spin-up by transfer of angular 
momentum by the accreting material |01 • 

Analysis of the recent ToO observation of 
HETE J1900. 1-2455 is also discussed. We observe 
in this object Type I X-ray bursts showing, at high 
energy, evidence of photospheric radius expansion in the 
burst profile. Using JEM-X/ISGRI data at high energy, 
and assuming the bolometric burst peak luminosity 
during photospheric radius expansion to be saturated at 
the Eddington limit, we measure the source distance to 
be around 5 kpc. A very interesting result is that the 
pulsations at the spin frequency are not always observed 
in this source. 



2. X-RAY SPECTRA 
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Figure 1. Simultaneous INTEGRAL, XMM-Newton and 
RXTE spectra ofXTE J 1807-294 fitted with an absorbed 
disk black body, DISKBB, plus COMPPS model. The 
EPIC-pn and M0S2 spectra in the 0.5-10 keV range and 
PCA/HEXTE in the 3-200 keV and IBIS/ISGRI in the 20- 
200 keV range are shown. The DISKBB model is shown 
by a dotted curve, the dashed curve gives the COMPPS 
model, while the total spectrum is shown by a solid curve. 
The lower panel presents the residuals. 

For the INTEGRAL observed MSP sources the high en- 
ergy spectrum could be described with a photon index of 
- 1.8-2.0 and a cut-off energy at - 80-100 keV. How- 
ever, the model does not describe the spectrum below 
15 keV well, which requires a more complex descrip- 
tion. MSP energy broad-band spectra are successfully 
fitted by a two-component model consisting of a multi- 
color blackbody soft X-ray emission and a Comptonized 
spectrum, for the hard X-ray emission. The soft thermal 
component could be associated with the radiation from 
the accretion disc and/or the heated NS surface around 
the shock. The hard emission is likely to be produced 
by thermal Comptonization in the hot accretion shock on 
the NS surface with the seed photons coming from the 
stellar surface. The spectra are best fitted with the ther- 
mal Comptonization model COMPPS in the slab geome- 
try rt l9[l . The main model parameters are the Thomson 
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Figure 2. The unfolded spectrum of ICR J0029 1+5934 
fitted with an absorbed COMPPS model. The data points 
correspond to the JEM-X (5-20 keV) and ISGRI (20- 
200 keV) spectra, respectively. The total spectrum of 
the model is shown by a solid curve. The lower panel 
presents the residuals between the data and the model. 

optical depth tt across the slab, the electron temperature 
fcTc, the soft seed photon temperature kTsccd, and the in- 
clination angle 6 between the slab normal and the line of 
sight. The seed photons are supposed to be injected from 
the bottom of the slab. The soft thermal emission is fitted 
by a simple blackbody, BBj or a multi-temperature disc 
blackbody, DBB, model [18]. The best fit parameters for 
the three observed MSP with INTEGRAL are reported in 
Table [T] In Fig. [Tl|2] and [3] we show the unfolded spec- 
trum and the residuals of the data to the BB or DBB plus 
COMPPS model. 

The spectra of XTE J 1807-294, IGR J00291H-5934, and 
HETE J 1900. 1-2455 are well described by a combination 
of thermal Comptonization and a disk black body or sim- 
ple black body. The hard spectral component contributes 
most of the observed flux (70-80 per cent), even though a 
soft component (disk black body) is required by the data. 
If the accretion disk is truncated at the magnetosphere ra- 
dius where the material can accrete along the magnetic 
field towards the pole of a NS, the matter over the pole is 
heated by a shock to a temperature of ^ 30-50 keV. The 
hot spot at the surface with temperature ^ 0.8-1.5 keV 
gives rise to the seed photons for Comptonization in the 
hot plasma. Since the hard X-ray emission is pulsed (see 
Sec. 3. 1), a fraction of it must originate from the regions 
confined by the magnetic field. The most obvious source 
of hard X-rays is the place where material collimated by 
the magnetic fields impacts the surface. 

For XTE J 1807-294, if we consider the inclination of the 
system to be 60° < i < 83°, this allows us to deter- 
mine the inner disk radius which lies in the range 20- 
40 km (for the distance of 8 kpc). The apparent area of 
the seed photons is Asccd ^ 26{D/8 kpc)^km^, which 
is also what expected from a hot spot on a NS surface. 
Many characteristics of the source are similar to those 
observed in other MSPs, SAX J1808.4-3658 id S 
and XTE J1751-305 lEl. For IGR J00291H-5934 no 
additional thermal-like component (either blackbody or 
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Figure 3. The unfolded spectrum of HETE J 1900. 1-245 5 
fitted with a COMPPS model plus a BB and a gaussian 
line. The data points correspond to the PCA (3-22 keV), 
JEM-X (3-22 keV), HEXTE (16-90 keV) and the ISGRI 
(20-300 keV) spectra, respectively. The blackbody model 
is shown by a dot-dashed curve, the dotted curve gives 
the COMPPS model, the dashed curve is the line, the total 
spectrum is shown by the solid curve. The lower panel 
presents the residuals between the data and the model. 



temperature depends on the optical depth, but tt x kT^ 
is approximately constant. 



3. TIMING 

The timing analysis of IGR J0029 1+5934 showed for 
the first time that INTEGRAL/ISGRI is capable of de- 
tecting the pulse profile of a 1.67 ms pulsar up to 150 
keV. This allowed us to study the pulsed fraction up to 
this high energy. We confirmed for the first time the 
increase of the pulsed fraction with energy in an accre- 
tion powered MSP. This can be explained by two mod- 
els: the action of the Doppler effect on the exponentially 
cutoff Comptonization spectrum from the hot spot, or a 
componization model where the pulsed fraction is pro- 
duced in a corona cloud. We also measured soft time 
lags with a complex energy dependence. Similarly to 
SAX J1808.4-3658 and XTE J1751-350, the time lags 
increase rapidly with energy until 10 keV. However, in 
IGR J0029 1+5934, the time lags first increase and then 
decrease slightly, saturating above 15 keV, and possibly 
reaching zero around 50 keV, see Fig. 6. 



disk blackbody) was required by the fit, very likely be- 
cause the bulk of its emission occurs below 3 keV, out- 
side the covered energy range. However, the apparent 
area of the seed photons, which turns out to be ^scod ~ 
2\{D /b kpc)^ km^, could correspond to a hot spot radius 
of 2.5 km during the INTEGRAL ToO outburst phase. 

Only HETE J1900.1-2455's spectrum differs from other 
MSP in requiring thermal soft X-ray emission with nearly 
double the temperature. We could not distinguish be- 
tween the multi-temperature blackbody and single black- 
body models, as both gave comparable parameters and 
. However, for a distance of 5 kpc, the disk blackbody 
gives an inner disk radius, i?in Vcos7 = 2.6 km, smaller 
than the expected NS radius. For the blackbody emis- 
sion, the fit implies the apparent area of the emission re- 
gion Asccd ^ km^, which could be consistent with the 
heated NS surface around the accretion shock |20, T6ll. 
From our spectral fits we infer that this emission is not 
likely produced in a multi-temperature accretion disk but 
more likely arises from thermal emission at the NS sur- 
face. 

We note that the product tt x kTc is very close in the three 
sources. The spectral shape and the product tt x kT^ 
are very stable during the outburst as is observed in other 
sources too [e.g. 16]. The constancy of the spectral slopes 
during the outbursts and their extreme similarity in differ- 
ent MSPs can be used as an argument that the emission 
region geometry does not depend on the accretion rate. If 
the energy dissipation takes place in a hot shock, while 
the cooling of the electrons (that emit X/7-rays via ther- 
mal Comptonization) is determined by the reprocessing 
of the hard X-ray radiation at the neutron star surface (so 
called two-phase model, the spectral slope is determined 
by the energy balance in the hot phase and is a function 
of the geometry. At constant geometry (e.g. slab), the 



3.1. Pulsed fraction of X-rays 
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Figure 4. This cartoon illustrates the different emission 
patterns responsible for the time lags of the pulsed emis- 
sion. 0g°* and are the dimensionless temperature 
of the accretion column and reflector, respectively, and 
Ciw ^0.1 is the disk illumination fraction. Soft time lag 
of the pulsed emission is the result of downscattering of 
hard X-ray photons in the relative cold plasma of the disk. 
A fraction of hard X-ray photons 1 — Cm directly seen by 
the observer are upscattered soft photons coming from 
NS and the disk. We account for the pulsed flux (pulsed 
fraction) as being produced in the corona cloud. 

We found for IGR J0029 1+5934 using the RXTE and IN- 
TEGRAL/ISGRI data that the pulsed fraction gradually 
increases with energy from ~ 6% at 6 keV to ~ 12 — 20% 
near 100 keV. This is the first time such behaviour has 
been measured for any of the accretion-powered MSPs. 

Since MSPs are rapidly rotating, a first interpretation is 
given through a Doppler boosting model which affects 
the fast rotating spot emitting pattern as a black body. 



Table 1. Best-fit spectral parameters with COMPPS + BB for DISKBBj model. 
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When the spot moves towards the observer, the emission 
increases, while for a spot moving away, the flux drops. 
The Doppler factor reaches the maximum a quarter of a 
period before the peak of the projected area, shifting the 
emission peak towards an earlier phase. The observed 
flux due to the Doppler effect varies as the Doppler factor 
to the power (3 + F) |20, 26], where the photon index F 
could be a function of energy. If the Doppler factor varies 
around 1 with 2% amplitude, we get 10% variability. 

The Comptonized spectrum can be approximated as 

Fe oc S-(ro-i) exp {^[E/E,]^) , (1) 

where Ec ^ kTc is the energy of the cutoff and parameter 
(3^2 describes its sharpness. The local photon index is 
then 

r{E) ^ 1 - = Fo + f3{E/E,f. (2) 

d m ii 

At low energies, F w Fq, and rms (or pulsed fraction) 
is a very weak function of energy. Close to the cutoff, 
the spectral index rapidly increases and the pulsed frac- 
tion should grows with energy, as observed (see Fig.|5t). 
The Comptonization models predict softening of the to- 
tal spectrum with simultaneous hardening of the pulsed 
spectrum at higher energies. 

On the other hand, the pulsed fraction of 
IGR J0029 1+5934 vs energy can be also explained 
by the energy dependent electron cross-section 
(Te{E) = (Tt(1 — 2z) and consequently by Comp- 
ton cloud optical depth as a function energy 



Tel(^) =TT,cl(l-22), 



(3) 



where z = E/rrieC^ is a dimensionless photon energy. 
In Figure |5] we show the energy dependent pulsed frac- 
tion of (a) IGR J0029 1+5934, (b) XTE J 175 1-305 and 
(c) SAX J1808.4-3658. In fact, we assume that the accre- 
tion column is embedded in Compton cloud as shown in 
Figure |4] [see also 24, for the geometry details]. Because 
Tci {E) decreases with energy, at higher energies a larger 
fraction of the pulsed direct hard X-ray radiation 

Amsxs(-B) = Arms,os(0) exp[-rci(£;)] (4) 
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Figure 5. The observed energy dependent pulsed frac- 
tion of (a) IGR J0029 1+5934 (b) XTE J175F305, and 
(c) SAX J1808. 4-3658 along with the best-fit comtoniza- 
tion model. 



originating in the accretion column can escape to the ob- 
server. 

A different scenario of the energy dependent amplitude 
formation would be if one can assume that there is no 
electron (Compton) cloud between the accretion column 
(where the Comptonization spectrum is formed) and the 
observer, but the energy dependence of the amplitude 
is formed as a result of upscattering. The soft photons 
upscattered off of hot electrons (of the accretion col- 
umn) increase their energy with a number of scatterings. 
On the other hand, the amplitude of the pulsed radia- 
tion exponentially decreases with a number of scatterings 
^rms(i*) = ^rms(O) exp(— where (3 is the inverse of 
an average number of scatterings [see e.g.|23t] and conse- 



quently with a energy of the upscattered photon 



A. 



rms,up [ 



rms.up 



(5) 



(where ad — (3/49 is approximately the energy index 
of the Comptonization spectrum). A characteristic seed 
(disk and NS) photon energy is £^<; < 5 keV and con- 
sequently, all photons at higher energies (> 5 keV) are 
produced by upscattering. The amplitude of upscattered 
photons (in which the energy is higher 5 keV) should de- 
crease with energy. 

In the general case, some fraction of the pulsed soft pho- 
tons, A 



rms.up 



upscatter off hot electrons in the accretion 
column on the way out as another fraction of the pulsed 
photons, Amis. OS, forming the hard X-ray tail escape to 
the observer passing through the Compton cloud. In this 
case a following combination of Eqs. (|4]and|5]l 

A„„3(£^) = A,„3,up(£^/£^s)""='+A,„3,cse-"='(^' (6) 
leads us to the formula of the emergent pulsed amplitude. 

In Figure|5]we present the best-fit model (see Eq. |6]l along 
with data points. We fit only the data points which cor- 
respond to energies higher than the seed photon energy. 
We found for (a) IGR J00291H-5934 the Compton cloud 
optical depth to be rT,ci = 3 ± 0.2 and Arms.os = 98lg 
%, Ams.up = 5.2 ± 0.9 %. For (b) tt,c1 = 0.25 ± 0.09 
and Ams,cs = 3.5 ± 0.8 %, Ams,up = 1-8 ± 1.1 % 
and (c) rT,ci = 0.4 ± 0.2 and A„„s,cs = 5.9 ± 6 %, 
Ams.up = 2.5 ± 1.2 %, respectively. For all of the fits 
ad was ~ 0.8. 



3.2. Time lag 

Additional information for the X-ray production pro- 
cesses and emission environment can be obtained by 
studying the pulse profile and phase shift between X- 
ray pulses at different energy ranges. For SAX J1808.4- 
3658 and later for XTE J1751-305 it was found that the 
low-energy pulses lag behind the high-energy pulses (soft 
phase/time lags) monotonically increasing with energy 
and saturating at about 10-20 keV lH [Toj [l6ll . 

First, this phase/time lag effect was interpreted as a result 
of the photon delay due to downscattering of hard X-ray 
photons in the relatively cold plasma of the disk or NS 
surface 0, 0] . They argued that the photon time lags 
were an intrinsic signature of interaction of the Comp- 
tonized radiation with the NS and accretion disk plasma. 
Moreover, the absolute values of time lags (about hun- 
dreds /is) are consistent with the electron scattering time 
scale tc = tt{L/c). The effective Thomson optical 
depth of the cold reflector tt — nf-cr^L is about a few 
and typical sizes of the NS photosphere and half-width of 
the disk L are of the order of 10® cm. On the other hand, 
Poutanen & Gierlinski |20] suggested that the lags may 
be produced by a combination of different angular dis- 
tribution of the radiation components and relativistic ef- 
fects. The observed spectrum of IGR J00291H-5934 con- 
sists of a black body from the neutron star surface and 
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Figure 6. The measured soft time lag of the pulse pro- 
file versus energy (crosses) with respect to the first en- 
ergy channel. The best-fit curve using the Comptoniza- 
tion model (see Eq. 7} is shown with the solid line. The 
dot-dashed lines in panel (a) correspond to the upper and 
lower limits of the electron number densities of the Comp- 
tonization emission area, 71^°* and of the reflector, ri^^^ in 
IGR J00291+5934. The panels (b) and (c) are related to 
XTE J1751-305 and (c)for SAX Jl 808.4-3658. The dot- 
dashed lines correspond to the upper and lower limits of 



a component produced by Comptonization of these seed 
photons in the hot electron region, presumably a shock 
which can be represented as a plane-parallel slab. The 
angular distributions of the black body and Comptonized 
photons emitted by the slab are significantly different. 
The difference in the emission patterns causes the two 
components to show a different variability pattern as a 
function of the pulsar phase, with the hard Comptonized 
component leading the soft black body component. How- 
ever, the soft lags found up to ^ 100 keV in the 1.67 ms 
accreting MSP IGR J00291H-5934, show a more complex 
energy dependence [6]. 

Recently, Falanga & Titarchuck derived a Comp- 
tonization model for the observed properties of the en- 
ergy dependent soft/hard time lags. This accounts for the 
soft lag by downscattering of hard X-ray photons in the 
relatively cold plasma of the disk or neutron star surface. 



A fraction of soft X-ray photons coming from the disk 
or neutron star surface are upscattered off hot electrons 
of accretion column and this effect leads to the hard lags 
as a result of thermal Comptonization of the soft pho- 
tons. The model allows the observed soft and hard lags 
to be reproduced by the downscattered and upscattered 
radiation as a function of the electron number densities 
of the reflector, nl°^, and the accretion column, n^"*. 
In the case of the accretion-powered millisecond pulsars 
IGR J00291H-5934, XTE J1751-305, and SAX J1808.4- 
3658 the observed time lags agree very strongly with the 
model (see Fig. |6ll. The resulting time lags are seen 
as a linear combination of the positive (hard) time lags 
formed in the Comptonization emission area and the neg- 
ative (soft) ones formed in accretion disk and NS photo- 
sphere as a result of reflection of the hard radiation (see 
Fig.H: 



1 1 - A9,,t/z nf 1 - Ciii ln(z/z,) 
.4^?,ef ''l-40„f/2* nhot Cm ln[l + (3 + a)0hot] 

(7) 

where is the electron number density of the reflec- 
tor, ng°* is the electron number density of the Comp- 
tonization emission area (accretion column) and fljcf — 
kTl°^ / [meC^) is a dimensionless temperature of the re- 
flector. We assume a typical value of 6'iof < 0.7 keV/51 1 
keV. kTl^°^ and a are the best-fit parameters for the 
hot plasma temperature and spectral index of the Comp- 
tonization spectrum, and ^hot — kT^°^ / {nieC^). 

We also assume that the seed photon energy is near the 
lowest energy of the downscattered photons, i.e. Zg « 
z*. An observational value of E^, is about 3 keV, and 
we assume that the value of the illumination factor Cin 
is about 0.1 [see for details 9]. To examine the pre- 
dicted time lag using the Comptonization model Eq. (|7]), 
we use the observed pulse phase lag data and best-fit 
spectral parameters kT^°^ and a of the accreting MSPs 
IGR J00291H-5934, XTE J1751-305, and SAX J1808.4- 
3658. The SAX J1808.4-3658 data were taken from Cui 
et al. |l5|]. Ford jlOtl and represented in Figure |6] by 
open circles and filled squares, respectively. The data 
for XTE J 175 1-305 and IGR J00291H-5934 were taken 
by Gierliriski & Poutanen fl^, Falanga et al. |@]. For 
the reflector temperature we used fcT™' = 0.4 keV for 
all fits. The model to fit the time lags in XTE J1751- 
305 and SAX J1808.4-3658 has only one free parameter 
(see i^), the number density of the "cold" reflector n^g^. 
Presumably, in XTE J1751-305 and SAX J1808.4-3658 
the density of the Comptonization region is much higher 
than that of the "cold" reflector, i.e. n^°^ > n'-^°^. It 
is not surprising that the seed photons are Comptonized 
(upscattered) in the very dense plasma of the accretion 
column [see e.g. |2, [H. The best-fit values are n^"^ = 
6.3 X 10^^ cm-3 and 3.2 x 10^^ cm^^ for XTE J1751- 
305 and SAX J1808.4-3658, respectively. The fits of the 
time lag data for IGR J00291H-5934 provide us the best- 
fit values of the "cold" and hot plasma densities. Both 
the positive (upscattering) and negative (downscattering) 



time lags contribute to the apparent time lags because 
n'f = 6.9 X lO^** cm-3 and = 2.1 x lO^* cm^^ 
are of the same order of magnitude. 

The time lag data of IGR J00291H-5934, XTE J 175 1-305, 
and SAX J1808.4-3658 are collected in different time in- 
tervals lasting from hours to days [5, 10, 16, 6]. However, 
a hydrodynamical (density perturbation) time, thydro, in 
the innermost part of X-ray NS source is on the order 
of the ratio of the NS radius to the sound speed, namely 
thyd.ro ^ -RNs/csound ~ 0.1 s. Thus, during the data col- 
lection periods the densities of the surrounding plasma 
can vary. As a result, the time lags also change because 
they are very sensitive to density variations (see Eq. |7|. 
In fact, these density variations can be inferred from the 
time lag data. One can make a band between two curves 
of time lags vs energy which contain all of the time lag 
data points (see Fig. |6]). This allows constraints to be put 
on rt™' and n^°^ . The lower and upper curves in Figure 
|6]coiTespond (a) to nj^^^ = (6.1 - 8.0) x 10^^ cm"^ and 
nj;°* = (1.6 — 2.6) X 10^* cm^-^, respectively. For the 
sources (b) XTE J1751-305 and (c) SAX J1808.4-3658 
the density variations are n^"^ — (6.0 — 6.6) x 10^^ cm~'^ 
and nl°^ = (2.9 - 3.6) x 10^^ cm^^. Thus, the plasma 
density of the "cold" reflector can change as much as 10 
% during the entire data collection. 

Soft lags are observed only if n]i'^f <C n^"* and thus the 
relative fraction of the dowscattering time lags in the to- 
tal time lag sum is ^--^100%. However, Comptonization in 
a non-uniform accretion medium might account for the 
observed time lag as a non-monotonic function of en- 
ergy. Using this model, upper and lower limits of the 
atmosphere density (density variation) in the region of 
phase/time lag measurements were determined. Using 
the observed energy dependent pulsed amplitude one can 
find the variation of the Thomson optical depth of Comp- 
ton cloud in which the accretion column is embedded. 



4. BURSTS 



MSP are also known to exhibit Type-I X-ray bursts. 
These bursts were observed for SAX J1808.4-3658 and 
XTE J18 14-338 llllll. Several X-ray bursts have also 
been observed for HETE J 1900. 1-2455 by various obser- 
vatories 1 12]. One of these bursts was also observed dur- 
ing the INTEGRAL ToO observation [M- In FigureQwe 
show the JEM-X and ISGRI burst light curve (28 Octo- 
ber 2005, 10:25:12 UTC) in different energy bands. The 
burst rise time was 0.23 ± 0.05 s. The double peak pro- 
file is clearly evident at high energy (lower panel) within 
the first 12 s, while during this time the intensity at lower 
energy (upper panel) remains constant. This can be inter- 
preted as a consequence of a photospheric radius expan- 
sion (PRE) episode during the first part of the outburst 
[see e.g.. I25I1 . When a burst undergoes a PRE episode, 
the luminosity remains nearly constant at the Edding- 
ton value, the atmosphere expands, and its temperature 
decreases resulting in a double-peak profile observed at 



high energies. The tail of the burst at high energy can be 
seen for about 5 s after the PRE episode. 

A time-resolved analysis of the net burst spectrum based 
on the JEM-X/ISGRI 3-50 keV energy band data was 
well fit by a photoelectrically absorbed blackbody. Dur- 
ing the first 12 s, the unabsorbed bolometric flux was al- 
most constant at Fpcak = 9.5(2) x 10~^ erg cm^^ s^^, 
while the blackbody temperature dropped in the middle, 
simultaneously with an increase by a factor of ~ 1.5 in 
blackbody radius. The observed temperature reached a 
peak at ~ 2.5 keV, and then gradually decreased. The 
softening of the emission towards the end of the decay 
phase is also indicated by the e-folding decay times of 
12.5 ± 0.5 s in the 3-6 keV to 4.3 ± 0.7 s in the 12- 
20 keV energy band. This behavior is typically observed 
during PRE X-ray bursts. 

The burst fluence was /b = 1.67(6) x 10~^erg cm^^, 
calculated by integrating the measured Fboi.bb over the 
burst duration of ~ 50 s. The effective burst dura- 
tion was T = /bZ-Fpoak = 18.2(8) s, and the ratio 
of the observed persistent flux to the net peak flux was 
7 — impels /-Ppcak = 0.021(1). The burst has the same 
spectral parameters as previous bursts for this source ob- 
served with HETE-2 or RXTE fljil . Assuming a helium 
burst at the Eddington limit and canonical NS parameters 
(1.4 solar mass and radius of 10 km), the source distance 
is estimated to be ^ 5 kpc. 

From the observed INTEGRAL burst properties and the 
mass accretion rates inferred from the persistent luminos- 
ity, the present theory predicts that this burst is pure he- 
lium burning. For helium flashes, the fuel burns rapidly, 
since there are no slow weak interactions, and the local 
Eddington limit is often exceeded. These conditions lead 
to PRE bursts with a duration, set mostly by the time 
it takes the heat to escape, of the order of 5-10 s, as 
observed. In the framework of the thermonuclear-flash 
models [e.g., |25t] the burst duration, r < 20 s, and the 
ratio of the observed persistent flux to the net peak flux 
7 w 0.02 indicate a hydrogen-poor burst which is in 
agreement with the companion star being a helium-rich 
brown dwarf llTll . 

Because there were no other bursts observed during 
the INTEGRAL observation, the burst recurrence time, 
Atrec, must be at least one day. We can compute the 
ratio of the total energy emitted in the persistent flux to 
that emitted in the burst 



fb 



'-Mr 



7 



0.021 
18.2 



Mrec > 100, 



which is consistent with pure helium bursts [see e.g. 12511 . 
Taking the burst total energy release £'(, = 5 x 10^^ 
erg (derived from the fluence fb) and He burning effi- 
ciency of eHc ~ 1.7 MeV/nucleon w 1.6 x 10^^ erg g~^, 
we estimated the amount of fuel burned during the burst 
Eb /chc ~ 3.1 X lO'^^ g. For the mean mass accretion rate 
of 2 per cent of the Eddington, a burst recurrence time of 
2.2 days is expected. 
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Figure 7. A bright X-ray burst detected from 
HETE J1900.1-2455. The JEM-X (3-20 keV; upper 
panel) and IBIS/ISGRI ( 18^0 keV; lower panel) net light 
curves are shown (background subtracted). The time bin 
is 0.5 s for both IBIS/ISGRI and JEM-X light curve. At 
high energy the burst shows strong evidence of photo- 
spheric radius expansion. 



5. CONCLUSIONS 



The science of accreting millisecond pulsars has been 
carried forward by the high-energy capabilities of INTE- 
GRAL in concert with the timing resolution capabilities 
of RXTE and the spectral resolution of XMM at low en- 
ergies. INTEGRAL observations have contributed to our 
ability to characterize the high energy spectrum of these 
sources, and to study their interesting time-lag and pulse 
fraction behavior over a wide energy range, up to ~150 
keV. 

The reason for the lack of coherent pulsations in the per- 
sistent emission from LMXBs has been a longstanding 
open question in X-ray astronomy. For the first time, 
we are able to observe a system (HETE J1900. 1-2455) 
known to contain a millisecond period pulsar, which 
has continued in outburst after its pulsations have disap- 
peared. The transition of HETE J 1900. 1-2455 from an 
X-ray millisecond pulsar to a persistent LMXB could in- 
dicate that there is a population of suppressed X-ray mil- 
lisecond pulsars among the non-pulsating LMXBs. De- 
tailed observations of this source at the epoch of pulsation 
supression can help to solving the long-standing issue of 



missing pulsations in persistent LMXB emission. 
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